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v^-. ABSTRACT 
0^ 

We add to the lore of spherical, stellar-system models a two-parameter family 
Q^i with an anisotropic velocity dispersion, and a central point mass ("black hole"). 

The ratio of the tangential to radial dispersions is constant-and constitutes the first 
parameter-while each decreases with radius as r~^/^. The second parameter is the ratio 
of the central point mass to the total mass. The Jeans equation is solved to give the 
^ density law in closed form: p oc (r/ro)~'^/[l + (r/ro)^~^]2, where tq is an arbitrary scale 

^ ■ factor. The two parameters enter the density law only through their combination 7. 

^3 ' of the velocities having a Keplerian run, but with a variable anisotropy ratio. This 

o , 

\^ , law and the dispersion runs, which contain more than one radius-scale parameter. 
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O \ ellipticals-galaxies: kinematics and dynamics-galaxies: structure 



1. INTRODUCTION 

Model spherical stellar systems are useful in various contexts, as models of elliptical galaxies, 
star clusters, cores of galaxies, etc., and also as test beds for various theoretical ideas. To the 
library of classic models such as the Emden polytropes, isothermal spheres. King, and Michie 
models [see e.g. Binney and Tremaine (1987) and references therein], were added, in recent years. 



the models of Jaffe (1983) , of Bertin and Stiavelli (1989) , and of [Hernquist (1990) and their 



generalizations to "eta" models by [Dehnen (1993)| , and by p^remaine et al. (1994) 



There are various routes to constructing such spherical models: One may start from a 
distribution function (hereafter, DF) satisfying the Jeans theorem, and see if one gets a useful 
density and dispersion profiles. One is then ensured from the start that the model is an equilibrium 
one (but stability may still be a question). Trial and error is then required to obtain useful models, 
with relevant density and velocity-dispersion profiles. Alternatively, we may start by dictating 
constraints on the more directly observable density, p(r), and the runs of tangential and radial 
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velocity dispersions, crt{r) and ar{r), respectively. A necessary (but not sufficient) condition that 
there be an underlying stationary DF is that these three functions satisfy the Jeans equation 

1 2^ _ _M(r)G 

p ar r 

with ^ 

M{r) =4tt [ r'2 p(r') dr', (2) 
Jo 

and 

r2- 



/3^1-4- (3) 



This is still only one functional constraint on the three functions, so to pinpoint a model we need 
two further functional constraint. For example, we can restrict ourselves to isotropic models 
and put (7r{r) = crt(r) = cr{r). Then we can start with a density profile and calculate cr(r) from 
the Jeans equation, as in [Tremaine et al. (1994)] . We can, alternatively, assume an equation of 



state-a relation between p and cr-and solve the Jeans equation for both p{r) and cr(r) (as for gas 
polytropes). Or, we may start with an assumed o"(r) and solve the Jeans equation for p(r) (as in 
the case of isothermal spheres). 

If we want anisotropic models we may start with a different constraint on crr{r), crt{r). For 
example, we may assume that their ratio is position independent (as we do in much of this paper), 
and then dictate either a density or a dispersion profile. 

Such models are useful if at least some aspects of them can be given in simple, closed forms. 
In the end one has to ascertain that the model is underpinned by a positive DF satisfying the 
Jeans theorem, if one seeks to insure that the model is a realization of at least one equilibrium 
system. 

We offer here what we believe is a novel family of such models. We start by dictating the run 
of velocity dispersions: a constant ratio between the tangential and radial dispersions, while each 
of them decreases, in a Keplerian manner, as r^^/^. It is clear that in order to obtain a bound 
system with a finite total mass the dispersions must decline at least that fast asymptotically 
(lest they outstrip the escape velocity, which decreases in this manner). In fact, for a fixed-/? the 
dispersions must behave so asymptotically. Such a behavior is also "natural" for a stellar system 
near a central, dominant point mass. We then extend this family of models by assuming that only 
(cj^ + 0"^)^/^ is exactly Keplerian, but the anisotropy ratio is not constant. 

The useful features of the models are: simple, closed-form expressions for the density, 
accumulated mass, velocity dispersions, and potential, and the effortless inclusion of a central 
point mass. 

We describe the constant-/? models in section 2, derive some additional properties in section 
3 (including a partial discussion of underlying DFs), and explore a generalization in section 4. 
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2. THE MODELS 



We substitute our assumed dispersion profile, 



(l-/3)-V? = a2 = 5/r, (4) 



in the Jeans equation([l|) to obtain 

din p 
din r 



+ 2/? - 1 = -m(r), (5) 



where, m{r) = M{r)G/S is the dimensionless accumulated mass (M = SG ^ has the dimensions 
of mass and can be used to normalize all the masses in the problem), and (5 is constant. 

The addition of a point mass Mq at the center is effected by subtracting from the right-hand 
side of eq.(^) the constant 

mo = Mo/M = MoG/S. (6) 



We then obtain for the Jeans equation 

din p 



+ 7 = -m(r), (7) 



din r 
where 

7 = 2/3 + mo-l, (8) 

and m(r) includes only the accumulated mass of the stellar "gas", excluding the point mass; we 
thus have m(0) = 0. Equation(|7|) can be solved analytically: We eliminate p using 

M , 

to get an equation for m{r) {m! is the r-derivative of m) 

rm" + (7 - 2)m' + mm' = 0. (10) 

The left-hand side is the derivative of rm' -|- (7 — 3)m -|- which must then be a constant. 

Since m(0) = 0, this constant must be so we get 

-| . 1 , 1 < 1 , , 

m (3 — 7 m) m = —. (11) 

2 r 

For 7 > 3, m' is negative, which is unphysical. For 7 < 3, the equation is straightforwardly 
integrated to yield the desired mass profile 

"^^=2(3-7)^^^, (12) 

where x = r/ro, and the scale radius, ro, is an arbitrary integration constant; it equals the 
half-mass radius. 
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This mass profile corresponds to the density profile 



(3 - jfM x-^ 



Near the origin, the density behaves as r ^ while at infinity it behaves as r' 
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We see from eq.(12) that the total dimensionless mass is m(oo) = 2(3 — 7), or 

GMg = 2(3 - 7)5 = 4(2 - ())S - 2MoG. (14) 

S is thus determined by the two masses: the total "gas" mass, Mg, and the central point mass, 
Mo, through 

S = [A{2-/3)]-^G{Mg + 2Mo). (15) 

This is the virial relation for the models. We shall see in section 3 that for 7 > 2 the kinetic and 
potential energies of the model are infinite, but relation(|l5|) still holds for all legitimate values of 
7 < 3, and follows directly from the structure equation. 

The two masses determine the parameter uiq through tjiq = 4(2 — (3)Mo/{Mg + 2Mq), and 7 
through 

7 = 2/3 - 1 + 4(2 - P) ^°/^f . (16) 
' ^ ^ ^'l + 2Mo/Mg ^ ' 

In the limit of test-particle "gas" , Mg — > 0, with (3 fixed, we have 7 ^ 3, so the model becomes 
meaningless. A meaningful test-particle model is obtained if we let /3 ^ —00 (circular orbits) such 
that 7 < 3; (3Mg is then fixed at (7 - 3)Mo/2. Then 5^0, but 2cr|r -2/35 GMq, and we 
get a model with test particles on Keplerian circular orbits around a point mass, the likes of which 
can be built with any density distribution. If we admix test particles on non-circular orbits the 
Keplerian fall-off which we require cannot be maintained. 

In the limit /3 — > — 00, but Mg and Mq kept constant, we have 7 ^ — 00, and the sphere 
becomes an infinitely thin shell at an arbitrary radius tq, with the particles moving on circular 
orbits: A self gravitating sphere cannot consist of particles on circular orbits with Keplerian 
velocities unless they are all at the same radius. 



For 7 = 2, one obtains the density distribution of the Jaffe model | Jaffe (1983) ] [which is the 



same as the "eta" model of [Tremaine et al. (1994) , and Dehnen (1993) with rj = 1]. This value 



of 7 can be obtained only with a central mass as the maximum value of 7 for niQ = is 1. In the 
isotropic case it corresponds to M^/Mg = 3/2. 



3. GENERAL PROPERTIES 
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Energies 

We now consider tiie potential run of the model. The point mass makes the usual contribution 
to the potential; that of the stellar "gas", ^pg, can be straightforwardly integrated from the 
expression 

d^g _ GM{r) _ GMg x^-T 
dr Tq 1 + x^~'^ 

to give 

Vair) = 2Fi[l, (3 - 7)"^ (4 - 7)(3 - 7)"^ -{r^/r)^'-''\ (18) 

where 2F1 is the hypergeometric function [see e.g. IGradshteyn and Ryzhik (1980) |. As 

GWl 

2-Fi(a, 6; c; 0) = 1 we have asymptotically ipg{r) — > as expected. 

As to the behavior of (pg at the origin, we use the behavior of 2F1 at large values of its 
argument to deduce that for 7 < 2 

GMg^ ^2-7\ ^ ^4-7 



For 7 = 2, diverges logarithmically at the origin, and for 2 < 7 < 3, yjg diverges as r'^~'^ . 
Some special cases: For 7=1, one needs 2-^1(1, 1/2; 3/2; — x^^) = x tg~^(l/x) to get 

^,(,) = -^<.-'f?). (20) 



For 7 = 2, 2-^1(1, 1; 2; —x~^) = x ln{l + 1/x), yielding the known expression for the Jaffe model 
^g{r) = {GMg/ro)ln[r/{r + tq)]. 

The total kinetic energy of a model sphere is infinite for 7 > 2, and for 7 < 2 is 

SMg3-2/3^ 



(1 


-7^ 




"1 








-7/ 



Ek = ^^r(^]r(^]. (21) 



The potential energy of the stellar "gas" is, for 7 < 2, 



£^ = _ f^) r fi^) (^ + (22) 



2ro V3-7/ V3-7/\3-7 Mg I ' 



By eq.(15) they satisfy E'fc = —Ep/2. 
Projected properties 
If we define the integrals 

noo ^-ti(^2 _ \2\-l/2 

/(/^,-;A)^y^ ^^^^^^^> dx, (23) 
then the projected surface density S(a), at a projected radius a = Arg, is given by 

S(a) = ^%(7-1,3-7;A). (24) 
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The line-of-sight (projected) velocity dispersion, cr\\{a), is 



"II ^""^ = V, /(7-l,3-7;A) • 



At large projected radii (a ^ ro) we can write 

3-7 

4^1/2 ^2 ^^^^ 



na) » 5^^4^A-^ (26) 



and the line-of-sight dispersion 

p2 f e-'Y 



(«) - ^ \i , ( 1 - -■ (27) 



As examples, for 7 = 1 we get asymptotically 



9 , , Qtt , 5 >S' , , 

32(1-5/')-. (28) 



for 7 = 2 

and for 7 PS 3 



2,, = 4,.,S 



<^f(«)« 3^(1 -5/3)-. (29) 



2 / N 71" 3 S* 



4(1-4/3)-, (30) 



Near the origin (a ^ ro) we can write, when 7 > 1, 



and 



For example, for 7 = 2 



3-7Mg r(^') 

aJ(a)«-(l-|/3)-. (33) 
" TT o a 

For 7 ~ 3, the behaviors near the origin and at large radii are the same. We see that unlike the 
density, and surface-density, run, the projected-dispersion run depends on both (3 and 7. For 

7 > 1, o"|| declines as r~^/^ at both ends; the proportionality constant is, in general, different, but 
may be the same for special choices of P and 7 (for example, for 7 = 2, /3 = 5/6). 

Underlying distribution functions 

We have not been able to ascertain that all the models discussed here have underlying 
non-negative DFs that satisfy the Jeans theorem. For the isotropic case, (3 = 0, one can try to 



-7- 



search for a DF that depends on the particle energy alone. The procedure is then straightforward. 
First, we note from eq. (4-139) in Binney and Tremaine (1987) that for this limited class of DFs the 
density must not increase with radius anywhere. This necessary condition excludes the (isotropic) 
models with 7 < 0, for which p increases near the origin. The DF is uniquely determined by the 
density distribution through the Eddington relation [eq.(4-140b) of Binney and Tremaine (1987) | 



/(e) 



1 



8i/2vr2 



d^T^ 



1 /dp 



(34) 



Here, e = —E, where E is the energy, ^ = — and p is viewed as a function of ^ as both are 
functions of r. The value ^' = occurs at radial infinity, where p oc r'^~^ decreases faster than 
^ oc r~^; thus, the second term in eq.(|3^ vanishes. If we change the integration variable to r, 
eq.(^4|) becomes 

1 d?p d"^ dr 



One can show that 



dV 

dr 



fie) 



m + 7 

) 

m + 7 + 1 



^1/2, 



7 — 1 + m 



+ m(27 + 3/2) + (7 - 1)(7 + 3) 
(m + 7)(m + 7 + 1) 



(35) 



(36) 



Thus, for 7 > 1 the integrand is non-negative, and the DF is always positive. For 7 < 1, the 
integrand becomes negative for small enough radii, and the DF becomes negative for high enough 
energies. Such models do not have legitimate DFs that depend only on E. 

The two limiting models discussed at the end of section 2 clearly also have legitimate 
underlying DFs. 



4. A GENERALIZATION 



Scott Tremaine (private communication) has pointed out to us that the crucial step in our 
derivation, of once integrating the Jeans equation [going from eq.(|lO|) to eq.([ll|)], does not require 
that (Tr and at are separately Keplerian, it is enough that (erf +crf)^''^ is. Making only this relaxed 
assumption permits us to explore a larger family of models, which we now proceed to do. 

The general Jeans equation(||), in the presence of a central point mass Mq, can be brought to 
the form 

{r^M'^y - \2r{al + al) - GM^M' + -(M^)' = 0, (37) 
where the derivation is with respect to r. Now assume only that 



2,2 c -1 
+ = Sr , 



(38) 
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and, as before, define the dimensionless quantities niQ = GMq/S, rh{r) = GM{r) / S , and also 



s{r) = al{r)T/S = if^y (39) 



Then, the Jeans equation becomes 

(rsm')' - 7?m' + ^(m^)' = 0, (40) 

where t] = 2 — rriQ. As explained in section 1, having imposed only one functional condition on 
the three functions p, Cr, and at, we have yet to specify another in order to pinpoint a model. 
This can be a dictation of the run of the anisotropy parameter, as we shall do below by dictating 
s{r). Note that < s{r) < 1, where s = 0, 1 for purely tangential, and purely radial dispersion, 
respectively. 

As before, the left-hand side of eq.(^) is a derivative, and is readily integrated once. As s is 
bounded, the integration constant vanishes, and we get 

rsrh' — rjm + 2"^^ ~ ^- (^^) 

Values of r/ < are non-physical as they give a negative m'. For t] > 0, the integration of eq.(^lj) 
gives the mass profile 

m(r) = ' ^, (42) 

where 

Xir)=r^, (43) 



ro 



rs(r) ' 



and ro is some arbitrary radius that will appear as a scale in the density law. The integration 
constant in eq.(^) is chosen so that rg is the half-mass radius. Since s cannot exceed unity, X{r) 
must diverge (at least logarithmically) for r ^ 00, and so rh(oo) = 2rj, yielding the virial relation 

4S = G{Mg + 2Mo). (44) 

Thus, r/ is determined by 

77 = 2 . (45) 

The results for the constant-/? case are reproduced when we note that in this case S and S are 
related by 5 = (2 - P)S, so 777 = 777/(2 - /?), 7? = (3 - 7)/(2 - /?), etc.. 



The resulting density run is 



^^"■^ ~ 27rr3s(r)(l + e'?^)2- ^ ^ 



At large r, where X(r) is positive and diverges, the asymptotic form of the density is 

p(r^oo)«-^^e-^^; (47) 
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near the origin, where X{r) also diverges, but is negative, 

One can try different forms of s{r), subject to < s{r) < 1, and see if any interesting mass 
distributions result. Some general points first: If we depart from the constant-/? models discussed 
above, a new radius scale must be introduced by the choice of s, beside the scale, ro, introduced 
through the integration constant. This is because s must be of the form s{r/a). 

Near the origin the behavior is like that of the constant-/? models with (3 = /3(0). Thus, if 
s(0) 7^ [/3(0) 7^ — oo], X{r) diverges logarithmically (and is negative) at the origin, and we get 
from eqs.(^3|)(|48|) a power-law density behavior: p oc r'^^^~''^/^^^'^\ If s(0) = 0, X(r) diverges faster 
than logarithmically, and p vanishes non-analytically there [as in the constant s(r) = case, where 
the "gas" is concentrated in an infinitely thin shell]. 

At large radii, if s(cxd) 7^ 0, X diverges logarithmically, and p has a power-law behavior: 
p oc r~[^^''/'**^°°)l . If s(oo) = 0, X diverges faster, and p vanishes faster than a power (e.g. as an 
exponential of a power). 

Of the many models that can be produced we consider one class in more detail: Take s(r) to 
vary monotonically between the values at r = 0, and r at infinity as 

s r = ; ' . 49) 

Here, d is a power that controls the sharpness of the transition of s from v to r, and, when 
either vanishes, d measures the speed with which s does {d —d is tantamount to r). We 
take v ^ Q {u = Q can be treated by changing the sign of d), then, for r 7^ 0, the integral in 
expression(p3|) for X gives 

x;^(i/ + rx|^)^/'^' ^ ^ 

where x = r/a, xq = r^/a, a = ri/v, and 5 = rj/r — rj/v. The mass profile is thus 

2r]x°'(u + Tx'^Y/'^ 

X^iu + TX^y/'^ + X'^il^ + TX'^Y/'^' ^ ' 

and evinces the expected appearance of two distinct scale lengths. The density law is obtained by 



substituting expression (50) in eq. 



r]Mg x'^{v^Tx'^fl'^x%{y^TxiYI'^ 
~ 27rr3s(r) \x'^(y ^ rx^)^^ + x'^{v ^ rx'^)'5/rf]2 ' ^ > 

The density is a power law in radius at both ends with a power —(3 — 77/1^) at the origin, and 
— (3 -|- T//r) at infinity. 
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The case r = requires special treatment: here one finds 



\XoJ 



(53) 



giving rise to a density run that decreases exponentially at large r: 



p{r) oc r 



{2+a)-{a/d){r/a)'i 



(54) 



for r — ^ oo. 



We thus get a class of models that are more flexible as regards the density run, but they 
contain two scale lengths instead of one, and are certainly less wieldy. 

I thank Scott Tremaine for useful comments and suggestions. 
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